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ABSTEACT 

This report examines the data from a proton spectrometer 
flown aboard Skylab The instrument is sensitive to protons 
in the energy range 18 to 400 MeV A partial failure at 
the start of the mission was traced to probable failure of 
the optical coupling between the Csl crystal and the photo- 
multiplier The loss of all proton counts after several weeks 
was due to the same cause combined with low temperatures 

The partial failure restricted spectral analysis to two 
energy bands, 18 to 27 MeV and 27 to 400 MeV The directional 
data showed that a Gaussian angular distribution parameter of 
at least 70*^ is required for the low energy band and at least 
40^ Lor the high energy band The data, integrated over angle, 
indicate that the APS model extrapolated down to 18-27 MeV 
IS high by factors of 2 to 5 over most of the B-L space mapped 
In the 27 to 400 MeV range, the APS model is 20 to 100 percent 
low at low and high values of L, and is high at medium L 
values in the B-L space mapped 
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1 0 INTRODUCTION 


A proton spectrometer, sponsored by the Space Sciences 
Laboratory, George C Marshall Space Flight Center, NASA, was 
placed aboard the Skylab vehicle The Principal Investigator 
was Dr Godehart Guenther, University of Alabama in Huntsville 
This report presents an analysis of the data obtained during 
the iirst dozen days of the mission 

Proton radiation fluxes were sampled four or five times 
per day as the Skylab cluster penetrated the South Atlantic 
anomaly The orbit altitude was 450 km at 50 degree 
inclination 

The proton spectrometer has been described in other 
articles ^ ^ A brief description is given here for 

Camiliarization purposes 

The spectrometer is designed to discriminate against 
a large electron background A collimator and detector 
telescope arrangement limit the acceptance cone to a 22 5 
degree half angle as shown in Figure 1 Protons enter the 
instrument through an outer shield, SI, which reduces the 
low energy electron flux D1 and D2 are solid state detectors 
with a 10 nanosecond coincidence requirement These detectors 
define the acceptance cone and help ensure that single 
particles are being examined In addition, an energy thresh- 
old requirement of 400 KeV is imposed on D2 to discriminate 
against electrons 

A proton which penetrates D1 and stops in D2 has 18 MeV 
This IS the low energy boundary on the first proton channel 
A proton which penetrates a second shield, S2 , with less than 
250 KeV possesses 27 MeV This defines the upper energy 
boundary of the first proton channel 
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D1 - D2 - 03 - NOT D4 COINCIDENCE 


SAI 05G1 


Figure 1 Spectrometer Schematic 



Protons with energy greater than 27 MeV will deposit 
more than the 250 KeV threshold energy in the thallium- 
activated, cesium iodide detector, D3 Another coincidence 
between D3 and the two solid state detectors is required 
Protons between 27 and approximately 150 MeV will be stopped 
in D3 These protons deposit more than one MeV in 02, thus 
setting a "high" flag in the logic circuits Protons between 
150 and 400 MeV deposit 0 4 to 1 0 MeV in D2 and set a "low" 
flag in the logic circuits 

Light generated in the Csl D3 detector is transmitted 
Lo a photomultiplier tube through a light coupling material 
The amplified signal is sent to an analog-to- digital converter 
(ADC) which generates a wave train with the number of pulses 
proportional to the energy deposited Another circuit counts 
the pulses, examines the "high" "low" flag, and increments 
the counter in one of seven energy channels (channels 2-8) 

Spurious counts and energies entering from the side 
are discriminated against by requiring an anticoincidence 
Irom a plastic guard detector, D4 

Details of other special features such as an accidental 

f 1 2^ 

coincidence counter are available in previous reports ' ’ ' 
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2 0 FLIGHT HISTORY 


The proton spectrometer suffered a partial failure 
either during the launch or during the low temperature period 
aftet Skylab achieved orbit and before the instrument was 
turned on The meteoroid/ insulation shield around the Work- 
shop ripped off during the powered phase of the launch One 
ol the solar panel wings was also lost The other solar 
panel wing was jammed until the first Skylab crew managed 
to Iree it several weeks later 

Fortunately, some power was available from the Apollo 
Telescope Mount solar panels This power was rigorously 
conserved by shutting down all nonessential equipment including 
the proton spectrometer Because the spectrometer is mounted 
externally on the anti-solar side near the top of the Multiple 
Docking Adaptor (MDA) and because it is thermally isolated 
from the MDA structure, the temperature inside the instrument 
Cell to about -40°F 

Eventually, two decisions were made which-brought the 
instrument temperature up into its operating range One 
decision caused the spacecraft to be rolled and tilted in 
an effort to shade the Workshop by the ATM solar panels The 
Workshop temperature was climbing to about 120°F at that time 
This movement brought the instrument into the sunlight The 
other decision which helped raise the temperature was to 
activate the instrument because it used only 5 watts with 
the heaters on The temperature rose to 35 - 40°F Several 
weeks later after repairs were made and the spacecraft was 
placed in the solar-inertial mode again, the temperature fell 
below zero, suggesting that the heaters were undersized or 
that the insulation was damaged in a mapner not apparent to 
the crew during their fly^-around survey 

JI/— 
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REPRODUCIBILITY OF THE 
UKiljilNiUj riiOL lo rOOIv 


It was difficult to determine the condition of the 
spectrometer immediately after it was turned on Radiation 
was encountered during passage through the South Atlantic 
anomaly and the horns of the outer belt sporadically during 
the day These data were taped and telemetered to ground 
stataons at intervals The data were placed into a large 
data base where it could be accessed through the MOPS system 
Alter 24 hours the data were moved to permanent storage where 
access was difficult The spectrometer data were available 
in several displays and in several formats It was quickly 
determined that the D4 count rate from the guard scintillator 
and the D3 total dose from the cesium iodide crystal were 
behaving as expected Software errors prevented access to 
detailed data for several days After corrections to the MOPS 
software were made, a small quantity of data became available 
which indicated the high energy proton channel contained a 
suspiciously large number of counts relative to the other 
channels 

Several weeks after mission start, the spacecraft was 
placed in the solar inertial mode The temperature in the 
spectrometer fell to -20°F, and the proton counts ceased 
Other spectrometer functions such as D4 count and D3 total 
dose remained operational 

An attempt was made to analyze the early proton data 
Un 1 ortunately , the programs which had been designed to strip 
out the spectrometer data from the telemetry contained errors 
which rendered the formatted data meaningless The software 
group responsible for those programs had been reassigned to 
more urgent projects Eventually the programs were corrected 
By this time a torrent of data was coming in from Skylab, 
threatening to saturate computer facilities and personnel 
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A command decision was made to defer all reprocessing requests 
until after mission completion The data were reprocessed 
and made available for analysis early in 1975 

Initial examination of the data showed that it had not 
operated properly An informal fault analysis committee was 
convened by Dr Tom Parnell, Division Chief, Cosmic Ray Physics, 
Space Sciences Laboratory (SSL), MSFC Other members included 
Dr Godehart Guenther, Principal Investigator, John Watts, SSL, 
George Detko, SSL, and Charles Hill, Science Applications, Inc 

The fault was isolated to the optical coupling between 
the D3 ces3um iodide crystal and the photomultiplier tube 
This optical coupling is a transparent, rubber-like gel which 
' seals the two interfaces with the aid of a slight mechanical 
compression The compression was adjusted to give an operating 
range of 10 to 90°F It is probable that part of the seal was 
affected by the vibration of launch and/or by the cold soak 
suffered by the spectrometer 

Loss of part of the seal lowers the effective gain of 
the instrument by reducing the amount of light going to the 
photomultiplier tube The flaw is intermittent in nature, 
as evidenced by the behavior of the self~calibration mode 

I 

The self-calibration feature is based on an Americium 
source embedded in the D3 cesium iodide crystal Normally, 
the signal from this source is below the threshold set on the 
counting tram (ADC) During self-calibration , the gain is 
increased to such an extent that counts enter channels 2 and 3 
Counts are accumulated until either channel 2 or 3 overflows 
at 1024 counts If the other channel is below a predetermined 
value, the gain is stepped up or down towards the balance 
point The process is then repeated A total of 16 gain steps 
are available covering approximately 30 percent total gain 
change 
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During operation the gain changes would often he 
stable for a period of time, then change rapidly over a five 
minute interval The non-overflowing channel would fall from 
over 900 counts to very low values, usually about 100 but 
sometimes down to 5 The allowable gain change could not 
compensate for the intermittent loss of gain probably caused 
by defects in the optical interface Restoral of gain often 
took place within 5 minutes After several weeks, the count 
rate fell to zero 

An effort has been made to salvage usable data from the 
proton spectrometer Recalibration was not possible due to 
the fluctuating nature of the defect Eventually a two-band 
analysis was chosen because more detailed spectral analysis 
was impossible 

Thw low band energy range is 18 to 27 MeV Protons in 
this range may penetrate the first absorber and the D1 and D2 
solid state detectors which define the acceptance angle, 
but are stopped in the second absorber before reaching the 
D3 crystal The high energy band range is approximately 27 
to 400 MeV These protons stop in or completely penetrate 
the D3 crystal 

The low energy band is not significantly affected by the 
defect because the D3 crystal is not involved The low energy 
end of the high band is slightly affected The spectrometer 
IS designed so that 27 MeV protons entering the acceptance 
cone will deposit about 250 KeV in D3 A threshold discrim- 
inator IS set to reject signals below 250 KeV The gain loss 
during fluctuations is estimated to be about 50 percent Thus 
the lower limit may be approximately 27 25 rather than 27 0 
MeV This shift is comparable to fuzziness due to non- 
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due to angle of incidence The high energy end may be shifted 
somewhat more, perhaps 10 MeV This shift is comparable to 
other uncertainties such as the angle of incidence effect 
Also, the environment contains relatively few protons near 
400 MeV 

For the above reasons, the nominal energy limits of the 
upper band are unchanged Unfortunately, it is impossible 
to unscramble the 7 energy channels in the upper band so they 
are simply added together in this study 
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3 0 ANALYSIS 


The proton spectrometer is a directional instrument The 
proton flux is anisotropic with very small intensity parallel 
to the magnetic field lines and maximum intensity perpendicular 
to the field lines Therefore, the analysis requires the 
altitude, latitude, and longitude of the Skylab cluster at the 
time of a measurement, the pointing direction of the spectro- 
meter, and the direction of the field lines 

The experimental data is compared to the Vette APS proton 
f 3 ') 

environment modeP '^for omnidirectional flux above 30 MeV A 
computer program finds the APS model omnidirectional flux, 
applies a range of assumed pitch angle distributions to the 
fLux model, folds in instrument response, and computes the 
values which would be seen by an idealized proton spectrometer 
immersed in the APS model environment 

The direction cosines of the spectrometer pointing vector 
are (0,-l,0) in the MDA coordinate system This direction is 
related to a system rotating with the earth through the appli- 
cation of 6 rotation transformations 


(°\ 

F E D C B A I -ij 

Here (PSl, PS2, PS3) are the desired direction cosines of the 
proton spectrometer axis in the rotating earth coordinate 
system at the time of the measurement 

Transformation A relates the MDA system to the OWS 
system, accounting for bending at the flexible Airlock-OWS 
coupling The transformation is derived from Skybet tape Euler 
angles (El, E2, E3) The result of the first transformation 
IS as follows 
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PSOWSl = Cos El cos E2 sin E3. - Sin El Sin E2 
PSOWS2 = -Cos El cos E3 

PSOWS3 = -Cos El sin E2 sin E3 - sin El cos E2 

Matrix B transforms from the OWS system to the ECI (Earth 
Centered Inertial, of date 1950 0) system The Z-axis of the 
ECI system points along the north pole of the spin vector 
The X-axiS is toward the vernal equinox In terms of the 
SKYBET-furnished Euler angles, (Ul, U2, U3), the B matrix 
IS as follows Note that CUl = cos Ul and SUl = sin Ul, etc 

/ CU3 CU2 -CUl SU3 CU2+SU1 SU2 +SU1 SU3 CU2+CU1 SU2 ^ 


+SU3 


CUl CU3 


-SUl CU3 


' -CU3 SU2 CUl SU3 SU2+SU1 CU2 -SUl SU3 SU2+CU1 CU2 > 

Matrix C transforms a column vector in ECI (epoch 1950) 
to ECI (epoch 1973) coordinates It is a precession correction 
Matrix C is furnished by the Dudley Observatory, Albany, 

New York It was checked against Transformation 33, pp 
418-419, of Methods of Orbit Determination, Escobal , John Wiley 
& Sons, Inc , New York, 1965, and found identical to within 
two units in the seventh place after the decimal 


99998 

00514 

00223 


42912 

04016 

46372 


- 00514 
99998 

- 00000 


04016 

67880 

57433 


- 00223 

- 00000 
99999 


46372 

57437 

75032 


Matrix D corrects for nutation (wobble) from ECI (epoch 
1973) to ECT (Earth Centered True, epoch 1973) This small 
correction is based on observational data and may be computed 
only after the fact Matrix D is furnished by the Dudley 
Observatory 


99999 

00007 

00003 


99968 

37903 

,19983 


00007 

99999 

00001 


37907 

99972 


- 00003 

- 00001 


19975 

09138 
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Matrix E rotates the ECT coordinate’ system around the Z- 
axis so that the X-axix is moved from the vernal equinox to 
the Greenwich meridian 


/ cos 6 
E= l-sin 0 


sin 0 
cos 0 


0 = 99 6909833 + 36000 7689 TU + 00038708 TU 

TU = elapsed time in Julian centuries, or 

TU = (2441682 5 - 2415020 )/36525 

A half day has been subtracted to change from Greenwich noon, 
Jan 0, 1950, to 0 hr, 0 min, Jan 0, 1973 

Matrix P corrects from ECT (1973 0) to the time of 
measurement 


cos wt 
sin wt 


sin wt 
cos wt 


t = time since 0 hr , 0 min , Jan 0, 1973 
w = 26251 61453 radians/hr 

The SIX rotations shown above convert the PS axis from 
the local MDA system to the rotating earth system 

The counts in proton channels P1-P8 are screened for 
validity Total counts fewer than 5 distributed among the 
8 channels during a 1 2 second period are discarded Counts 
greater than 50 in a single channel indicate a probable error 
and that period is discarded The counts are normalized to 
counts per second and dead time corrections are applied 
These calculations are made at each measurement point, i e , 
every 1 2 seconds along the trajectory where sensible flux is 
encountered The measured values and calculated values are 
summed over time, then ratiod to provide spectral correction 


tactors 
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At the same time, the data are tabulated according to 
location (B, L) , pitch angle (a), energy, assumed pitch angle 
distribution (a), and orientation in the earth-fixed coordinate 
system With the aid of this data, detailed 

corrections to the APS model environment are derived "Best 
lit" energy-dependent pitch angle distributions are also 
obtained Finally, some information is derived concerning 
the east-west asymmetry, though this effect is blurred at 
higher energies due to acceptance in the backward cone of the 
PS above 100 - 150 MeV 

The Vette APS flux parameters, a and b, are found by 
table look-up using values of B and L from the SKYBET tape 

APS FLUX (>E) = aE"^ 

This form is differentiated, and the resultant is integrated 
over energy ranges corresponding to those of the eight proton 
channels The angle a between the local magnetic field line, 

B, and the axis of the PS primary cone is found using the 
scalar product The angle a is used in a table look-up to 
determine that fraction of the APS omnidirectional flux which 
would be observed by an ideal PS The measured angular 
response of the instrument is triangular, peaked on the axis 
and going to zero at ^ 22 5° The assumed pitch angle 
distribution is Gaussian 

-(90 -g)^ 

Q~2g^~ 

a V 2 tt 

The Idealized instrument reading is obtained by numerical 
integration and condensed into tables for this program Note 
that a determines the angular distribution of the flux A 
set of 10 values are chosen for a so that parallel calculations 
may be made 
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The flux values computed above will be denoted by 
Each quantity will be further identified by two sets of 
subscripts 

The first set of subscripts is described as follows 
’ ^ijklm 

where 

1 = B' (note prime) 

j = L' (note prime) 
k = a 
1 = E 

m = a 

Figure 2 shows that rectangular boxes in B - L coordinates 
would be wasteful of computer storage, or would have large 
variations in flux level within a box Therefore, a trans- 
formation is applied to obtain B'-L’ coordinates which are 
better suited here 



The B' box limits are 0- 005, 005- 01, 01- 02, and > 02 
The L' box limits go from 0 to 7 in steps of 1, then > 7 

The box limits in the a dimension are 0-30, then 10 
degree steps to 90 degrees Angles greater than 90° are 
reflected back to the first quadrant 


The box limits in the energy dimension are 18-27, 27-400, 
and 18-400 MeV 
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The calculated fluxes, have an extra dimension, pitch 
angle distribution a, and 10 parallel calculations are made 
for this quantity The assumed pitch angles parameters are 
5, 10, 15, 20, 25, 30, 40, 50, 60 and 70 degrees 

A second set of subscripts and tables is necessary for 
the east-west asymmetry study These subscripts are 

ijlnln2 ’ ijlmnln2 

where 


1 = B' 
3 = L' 

1 = E 


m 

nl 

n2 


a 


0 


LV 

LV 


The o subscript, k, is replaced by two subscripts, nl and n2, 
identifying equal solid angle boxes in the Local Vertical 
system Here nl labels the polar angle with boundaries at 
0, 60, 90, 120 and 180°, while n2 labels the azimuthal angle 
with limits at 330, 30, 90, 150, 210, and 270° The solid 
angle is subdivided into 4 x 6 = 24 equal size boxes Here, 
the measured fluxes, are tabulated in a 5-dimensional 

table 

2 

Values of $*, $*2, $ and f are determined for each 
measurement (1 2 seconds) and are accumulated in separate 
tables according to the two sets of subscripts described above 
Auxiliary tables are compiled to show the number of entries 
in each box of each table 


The tables described above are too lengthy to keep in 
fast storage Therefore, $+, the 10 $-'s, and the two sets of 
subscripts are placed in a buffer When the buffer is 
Cull, it IS dumped into mass storage After the input tape 
IS processed, the values are brought back into fast core and 
processed The subscripts, maximum value, and quantity 
represented are listed below 
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1 8 E 
J 8 L> 


m lOo 
n2 6 tf) 


LV 

1 4 B' nl 4 0^^ 

The measured and calculated data are sorted into tables 
as described above Later, the 8 energy channels are 
collapsed into the 3 energy bands described previously due to 
the partical spectrometer failure The cumulative tables are 
stored on tape in the following order 

Number of Entries Tables 

1 Number of entries per bin for the pitch-angle tables 
(8 E's, 8 L's, 4 B's, 7 a's = 1792 words) 

2 Number of entries per bin for the east -west tables 
(8 E's, 8L's, 4B's, 6 azimuthal angles, 

4 polar angles = 6144 words) 

Pitch-Angle Tables 

2 

3 Measured pitch-angle, $*, 1792 x 2 words ) 

4-13 Calculated pitch-angle for 10 assumed distributions, 

<l> , 10 X 1792 X 2 words (0 ,$ ^) 
ra m’ m ' 

East-West Tables 

14 Measured east-west table, 6144 x 2 words 

15-24 Calculated east-west tables for 10 assumed 

distributions, $ , 10 x 6144 x 2 words 

m 

The total table storage is 182528 words 

Values are accumulated over each pass through the South 
Atlantic anomaly, then added to previous accumulations 
Within each box, the summed values are 

T = S = 

T2 = 2$+^ S2 = 2f^ 
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The mean value of the measured flux within a box is 
T/n, where n is the number of entries in that box A standard 
dovaaLion is defined herein (incorrectly) to be identical 
with the root mean square deviation The reason for this 
simplification is to reduce the alternate computations required 
when n = 1 

(S D )^ = i 2($+--p-)2 
n ^ ^ 


= I - 5% 


Thus 


(S D ) - nT2-T /n and the fractional (S D ) 


as (F S D ) = '^nT2-T /T 

for the measured flux within a box Similar expressions 
may be derived for the model fluxes 


EEPEODTJCIBILrEy OF THE 
0EB5WAL PAGE IS POOR 
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p(E > 5o MeV) 
cm' - sec 
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Figure 2 


Flux Contour Map of the Sky lab Orbit in B - L Space 
235 n. m. - 50° , 
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4 RESULTS 


The results of the analysis are shown in Tables Al, A2 , 
and A3 Table Al gives the ratio, R, of the measured flux 
to the model flux within each B’, L’ , a(I,J,K) box Table A2 
shows similar ratios after weighted summing over the pitch 
angle, a Table A3 shows the ratio of measured to model flux 
coming from the east hemisphere and the west hemisphere within 
each B' , L' box 

Within each table, data are presented for 3 energy groups, 

18-27, 27-400, and 18-400 MeV The measured fluxes are 

(S') 

compared to the APS model fluxes^ '^which are intended for use 
above 30 MeV Therefore, attention should be focused upon the 
center band (27-400 MeV) An attempt was made to use the 
newer AP7 proton environment model^^^ (valid above 50 MeV), 
but this model does not extrapolate well down to 27 MeV It 
may be noted that the APS and AP7 spectra, integrated over 
the entire orbit, are nearly identical above 50 MeV for the 
Skylab trajectory 

The data are taken from 3280 samples of 1 2 seconds 
each, or about 66 minutes of active measurement This data 
represents less than 2 percent of the data sought Of the 224 
B' , L’ , a boxes, 81 contain samples, averaging 40 samples per 


Table Al contains 10 sets of ratios The first, labeled 
"M=l”, assumes a 5 degree parameter in the APS angular 
distribution model This parameter increments by 5 degrees 
to 30, then 10 degrees to 70 The first set uses a flux 
model sharply peaked near the perpendicular to the magnetic 
field line This model predicts counts only when a is near 90 
degrees, K=5, 6, or 7 The ratio, R, of measured to model 
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counts xs often large for the K=5 boxes, whxch xndxcates the 
model counts are too small The zero values of R indicate that 
the ratios are so large that certain program constraints are 
exceeded The values labeled "ESDI" are fractional standard 
deviations (fsd) for the model counts The values labeled 
"PFSDl” are fsd's for the measured counts within a box 

Consider the 27-400 MeV results for M=1 (a=5*^) The 
model flux angular distribution is obviously too narrow The 
ratio R IS nearly reasonable only for K=7 and 6 where the 
spectrometer sees the maximum flux The R values range from 
0 286 to 1 267 for K=7 For K=5 (20-30° from the X) the model 
counts become too small and the ratio too large (11 to 953) 

At M=5 (0=25°), the number of zero ratios is 11 compared 
to 33 at M=1 Most of the other ratios lie between 0 5 and 3 
A few are still large, particularly for K=2 (30° <a<40°), 
where R lies between 2 2 and 23 2 It is evident that the 
real flux is again more smeared than the model flux 

At M=6 (o=30°), the number of zero ratios is 6 and 7 
values of R are greater than 6 0 Most of the large or zero 
(very large) values of R are for observations within 30 degrees 
of the field line The model is still predicting fluxes much 
smaller than those observed 

At M=7 (a=40°), the number of zero ratios in the 27 to 
400 MeV range is one and only 3 values of R are larger than 
3 0 The zero ratio box contains 2 samples The ratios larger 
than 3 0 contain 5, 6, and 15 samples 

This group of data shows some degree of consistency 
between the experimental and model data The assumed model 
data is smeared enough (n=40°) to simulate experimental 
measurements 

Smaller values of a cause large values of R at small 
pitch angle as shown in Figures 3, 4, and 5 Figure 3 is a 
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plot of R, the measured flux divided by the model flux, versus 
a for 1=1 (0<B'< 005) Figures 4 and 5 are for 1=2 
( 005<B’< 01) and 1=3 ( 01<B'< 02) Eight plots are shown 
in the figure, one for each L' band labeled J=l-8 These I, J 
bands may be correlated with the dashed boxes of Figure 2, 
which was constructed with the AP7 model rather than the APS 
model used herein V/ithin each plot the curves are labeled 
by the index K A value of K equal to 1 identifies measure- 
ments taken when the spectrometer view axis is within 30 


degrees of the field line, either parallel or antiparallel K=2 

, up to K=7 where 


corresponds to 30*^ to 40*^ look angles, etc 


the look angle is within 10° of perpendicularity and maximum 
count rates should occur 

The K=7 curves increase monotonxcally with o because 
the model flux (denominator of R) is decreasing The K=6 
curves exhibit similar behavior for large a, but may show a 
slight :^ise at small a For K=1 to 5, the value of R is 
very large at small values of a because the model flux is small 
at these angles 

The curves of Figures 3, 4, and 5 show that o should be 
at least 40° to 50° in the 27 to 400 MeV range in order for 
the model flux to simulate the measurements at small angles 
to the field lines 

Figures 6-8 show similar curves for the low energy group, 
18 to 27 MeV The absolute values of the ratio, R, are less 
significant because the APS model is not valid down to 18 
MeV However, the relative magnitudes of the K=l, 2 versus 
K=7 curves show that a a value of 70° or larger is necessary 
to permit the model flux to simulate the measurements in the 
low energy range 

The angular distribution of the trapped protons should be 

I 

energy dependent at 235 nautical miles altitude, because of 
the scavenging effect of the atmosphere Higher energy 
protons must have a smaller value of a. in the SA anomaly as 
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they approach their mirror point than the lower energy protons 
possess 

It IS unfortunate that the loss of spectral information 
prevents a better study of pitch angle versus energy 

Within each B' , L' box, the ratio of experimental to 
model counts has been summed over pitch angle, using weights 
determined from model fsd's The results are presented in 
Figures 9-11 Figure 9 is for the lower B' values. Figure 10 
IS for the middle B' values, and Figure 11 is for the higher 
B' values observed The top graph in each figure pertains to 

the 18-27 MeV band, the bottom graph, to the 27-400 MeV band 
Each curve is labeled by its L' band 

These figures must be interpreted carefully because some 
patch angle bands are empty Dashed lines indicate that some 
ratios have exceeded program constraints 

The low energy curves of Figures 9-11 are systematically 
small except for the top 3 L' bands for lowest B' At large 
values of a, the AP3 model-extrapolated down to 18 MeV- 
generally predicts more counts than the experiment measures 
by a factor of 2 to 5 The top three curves of Figure 9 show 
good agreement between experiment and the model in this region 

The high energy curves of Figures 9-11 should be 
interpreted at values of cr between 40° and 60° The curves 
are probably most valid for the low portion of the 27 to 400 
MeV energy range where most of the flux exists The ratio 
of measured to model counts is large for L' bands 1, 7, and 8 
The ratio is smaller for L' bands 2, 3, and 4 It would appear 
that the model fluxes are too low near the top and bottom of 
Figure 2 

The east-west asymmetry data of Table A-3 are not amenable 
to analysis because of insufficient data 
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Ratio of Measured to Model Counts in the 
East and West Directions as a Function of 
B ' and L ' 
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